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In a recent paper M. Dzugutov H describes a molecular dynamics cooling simulation where he 
obtained a large monatomic dodecagonal quasicrystal from a melt. The structure was stabilized 
by a special potential || designed to prevent the nucleation of simple dense crystal structures. In 
this comment we will give evidence that the ground state structure for Dzugutov's potential is an 
ordinary bcc crystal. More detailed publications of the phase diagram of the Dzugutov potential 
are in preparation [^,^) . 

The results are obtained by molecular dynamics simulation. The constraint method || has been 
applied to simulate NVT or NPT ensembles with an extension that permits constant tempera- 
^ ■ ture or pressure gradients also. The equations of motion are integrated by a fourth-order Gear 

ON ' algorithm. The time increment was <5i=0.005 (all quantities are given in Lennard-Jones units). 

The initial structures have been obtained by melting and equilibrating crystals with 250, 500 or 
1024 atoms at T=4.0 and densities of p = 0.63, 0.55 and 0.857 , resp. for 10 5 time steps at constant 
volume. The 250 atom sample has been cooled at rates of -0.002/5/; and -0.001/ With the faster 
rate one obtains a tetrahedrally close packed (tcp) structure Q where all atoms are 12, 14 or 15 
Q\ \ fold coordinated and surrounded by Frank-Kasper polyhedra. The quasicrystal also belongs to this 

class. The bcc structure is tetrahedrally packed to, but the coordination polyhedron is not of the 
Frank-Kasper type. The radial density functions of all the structures are quite similar, although 
the first maximum differs in shape between bcc and tcp in a good quality sample. Cooling the 
500 atom system at -0.001 / <5t and the 1024-atom system at -0.00025 / also leads to bcc crystals. 
At faster rates the large system shows partial ordering, but the structure as a whole did not form 
an ordered crystal. Since 250 and 1024 are magic numbers for a bcc crystal one could assume 
q \ that these samples get locked into the bcc structure. This is certainly not the case due to several 

reasons: first, the nucleating structures contain vacancies, which change the required number of 
atoms to a non-magic one. Second, the samples are sometimes twinned or contain more than one 
domain. Third, the basic [100] and equivalent lattice planes are never parallel to the simulation 
box, and the nucleated solids often contain screw dislocation loops and are globally twisted, but 
locally the atoms are in a bcc environment. 

Cooling the fluid at high constant pressure always yields a bcc structure which sometimes evolves 
into a fee structure at low enough temperature. At P=5 a tcp structure is obtained at -0.002/5t 
and a bcc crystal at -0.001/<5i. If the fluid is compressed at constant temperature the results are 
similar. In the range from T=2.5 to 1.0 the bcc structure forms at a compressing rate of Q.l/St. 
At T=0.6 a twinned crystal was found with a few tcp atoms. 

It has been shown by density functional calculations || and energy minimization calculations j^j 
that the bcc structure is the lowest in energy at densities where the pressure is low. At increasing 
density the tcp tr-phase || (a low order approximant of the quasicrystal) becomes more stable, 
and at high density the fee structure is the most stable. The difference between bcc and tr-phase 
is rather small so that it is currently not clear which one is more stable at low and moderate 
pressure and T > 0. 

At a high cooling rate one obtains an glass. In an intermediate range at low pressure a tetrahe- 
drally close packed structure or a quasicrystal is favoured since its structure is close to the fluid 
structure (Fig. [j]). At slow cooling rates the system has enough time to nucleate into the stable 
bcc structure. With increasing system size it takes more and more time for the fluid to reorder 
into the bcc structure, therefore the tcp structure will be generated. 
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FIG. 1. Bond angle correlation function for nearest neighbour atoms The histogram gives the freqeuncy 
of angles as a function of the cosine. The liquid (full line) is similar to a tcp structure (dotted line), whereas 
the perfect bcc structure is completely different, but close to the nucleated structure (dash-dotted line). 
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